Statistics of temperature fluctuations in an electron system out of equilibrium 
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We study the statistics of the fluctuating electron temperature in a metallic island coupled to 
reservoirs via resistive contacts and driven out of equilibrium by either a temperature or voltage 
difference between the reservoirs. The fluctuations of temperature are well-defined provided that 
the energy relaxation rate inside the island exceeds the rate of energy exchange with the reservoirs. 
We quantify these fluctuations in the regime beyond the Gaussian approximation and elucidate their 
dependence on the nature of the electronic contacts. 



The temperature of a given system is well-defined in 
the case when the system is coupled to and in equilib- 
rium with a reservoir at that temperature. Out of equi- 
librium, the temperature is determined by a balance of 
the different heat currents from/to the system How- 
ever, this applies only to the average temperature: the 
heat currents fluctuate giving rise to temperature fluctu- 
ations. Although the equilibrium fluctuations have been 
discussed in textbooks Q, their existence was still de- 
bated around the turn of 1990's 

In this Letter we generalize the concept of tempera- 
ture fluctuations to the nonequilibrium case by quanti- 
fying their statistics in an exemplary system: a metal 
island coupled to two reservoirs (see Fig. [TJ. The island 
can be biased either by a voltage or temperature differ- 
ence between the reservoirs. In this case, the tempera- 
ture of the electrons is not necessarily well defined. The 
electron-electron scattering inside the island may how- 
ever provide an efficient relaxation mechanism to drive 
the energy distribution of the electrons towards a Fermi 
distibution with a well-defined, but fluctuating temper- 
ature [l|, IH ■ Here we assume this quasiequilibrium limit 
where the time scale r e _ e of internal relaxation is much 
smaller than the scale te related to the energy exchange 
with the reservoirs. 

In equilibrium, the only relevant parameters character- 
izing the temperature fluctuation statistics (TFS) are the 
average temperature T a , fixed by the reservoirs, and the 
heat capacity C = ir 2 k B T a /(35i) of the system. The lat- 
ter is inversely proportional to the effective level spacing 
Si on the island. In terms of these quantities, the prob- 
ability of the electrons being at temperature T e reads 
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FIG. 1: Setup and limit considered in this work: a conducting 
island is coupled to reservoirs via electrical contacts charac- 
terized by the transmission eigenvalues {T"}. The temper- 
ature fluctuates on a time scale te characteristic for the en- 
ergy transport through the junctions. We assume the limit 
T e -e <C te where the internal relaxation within the island 
is much faster than the energy exchange with reservoirs. In 
this limit both the temperature and its fluctuations are well 
defined. 



this distribution we can for example infer the variance, 
((AT/) 2 ) = k B T 2 /C. As we show below, the scale of 
the probability log, InP ~ T a /8i is the same for the 
nonequilibrium case while its dependence on (T e /T a ) is 
essentially different. 

To generalize the concept of temperature fluctuations 
to the nonequilibrium case we examine the probability 
that the temperature of the island measured within a 
time interval tq . . . tq+At and averaged over the interval, 
equals T e : 
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P cq (T e ) oc exp — — = exp 

(1) 

corresponding to Boltzmann distribution of the total 
energy of the island. The probability has a Gaussian 
form even for large deviations from T a , apart from that 
the probability naturally vanishes for T e < 0. From 
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The average (•) is over the nonequilibrium state of the 
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system. The latter is evaluated using an extension of the 
Keldysh technique where the fluctuations of charge 
and heat are associated with two counting fields, x an d 
£, respectively @, @, The technique allows one to 
evaluate the full statistics of current fluctuations both 
for charge and heat current Q in an arbitrary multi- 
terminal system. In terms of the fluctuating temperature 
and chemical potential of the island, Tj{t) and /U/(t), and 
the associated counting fields £/(t) and x/(i) t ne average 
in Eq. ([2|) is presented in the form 

P(T e ) cx J Vb(t)TW I (t)T>xi(t)T>iJii(t)dk 



x exp < —A + ik 



tq+At 



dt(Ti(t)-T e ) 



(3) 



Here A = A[£i(t), 7>(i), xi(t), m(t)] is the Keldysh ac- 
tion of the system. The counting fields £j(i) and xi(t) 
enter as Lagrange multipliers that ensure the conserva- 
tion of charge and energy [9(. 

The Keldysh action consists of two types of terms, A = 
J dt{Si{t) + S c (t)), with Si(t) = Qixi+Eiii describing 
the storage of charge and heat on the island and S c de- 
scribing the contacts to the reservoirs. Here Qi — C c /ii 
is the charge on the island, Ej = C(T I )T I /2 + C c (j, 2 I /2 
gives the total electron energy of the island and C c is 
the electrical capacitance of the island. For the electri- 
cal contacts, the action can be expressed in terms of the 
Keldysh Green's functions as (we set /i = e = &b = 1 
for intermediate results) 
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The sums run over the lead and channel indices a and 
n. All products are convolutions over the inner time vari- 
ables. The trace is taken over the Keldysh indices and the 
action is evaluated with equal outer times. This action 
is a functional of the Keldysh Green's functions G a and 
Gi of the reservoirs and the island, respectively. It also 
depends on the transmission eigenvalues {7 1 "}, character- 
izing each contact. The counting fields enter the action 
by the gauge transformation of Green's function [§[ 



G{t,t') = e -5(x/-^(*)9 t )r 3( 5 (i,t') e l(xi+^(t')a t ')r3. 



where the Keldysh Green's function reads 
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For quasiequilibrium /(e) = {exp[(e — fi)/T] + is a 
Fermi distribution. In what follows, we assume the fields 
£(i),T(i) to vary slowly at the time scale T _1 , in which 
case we can approximate — i£(t)dt <— > £(£) e - 



The saddle point of the total action at x = C = 
yields the balance equations for charge and energy. As- 
suming that the electrical contacts dominate the energy 
transport, we get 



dQi 
dt 



dEj 



Cc d m =^Trf 3 ^T« 



[G a , Gj 



4 + T«({G/,G Q }-2) 



(7a) 



= CdtTj 



\G a , Gi 



4 + T£({G I ,G a }-2) 
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The right-hand sides are sums of the charge and heat cur- 
rents, respectively, flowing through the contacts a 

The time scale for the charge transport is given by 
r c = C e /G, with G = J2 a En T n/(^)- This is typically 
much smaller than the corresponding time scale for heat 
transport, t e — Ch/G t \ u where G t h = ir 2 GT/3. We as- 
sume that the measurement takes place between these 
time scales, t c <C At <C te- In this limit the potential 
and its counting field fii and xi follow adiabatically the 
Ti(t) and £i(t) and there is no charge accumulation on 
the island. As a result, we can neglect the charge ca- 
pacitance C c concentrating on the zero-frequency limit 
of charge transport. 

To determine the probability, we evaluate the path 
integral in Eq. ([3]) in the saddle-point approximation. 
There are four saddle-point equations, 
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Here Mb(t) — I inside the measurement interval (to, tq + 
At) and zero otherwise. Equations (|5a|) express the 
chemical potential and charge counting field in terms 
of instant values of temperature Tj and energy count- 
ing field £/, hi = /i/(£/,T/) xi = Xi((i, t i)- Tne tn i rd 
and fourth equations give the evolution of these variables. 
It is crucial for our analysis that these equations are of 
Hamilton form, £/ and Tf being conjugate variables, the 
total connector action S c being an integral of motion. 
Boundary conditions at t ± oo correspond to most prob- 
able configuration T e — T a . This implies S c = at tra- 
jectories of interest. 

The zeros of S c in £/ — T/ plane are concentrated in 
two branches that cross at the equilibrium point £/ = 
Q. T e = T a (Fig. a). Branch B = 0) corresponds to 
the usual "classical" relaxation to the equilibrium point 
from either higher or lower temperatures. Branch A cor- 
responds to "anti-relaxation": the trajectories following 
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FIG. 2: Time line of a huge fluctuation. The measurement 
is made at t ~ to with result T e = 2.5T a . For t < to, the 
temperature follows the "anti-relaxation" branch A, whereas 
after the measurement, it relaxes as predicted by a "classical" 
equation, (a) Contour plot of S c , evaluated for the case of two 
equal tunnel junctions with T\ = T2 = and bias voltage V 
(T a — v / 3eV^/(27rfcs)). The zeroes of S c are given by crossing 
curves; the dashed line indicates the boundary of the domain 
where £ is defined (£ > — l/(fcsT)). The sign of S c in the dif- 
ferent regions is also indicated. The arrows give the direction 
of the flow for the saddle-point solutions {£f (t),T?(t)}. (b) 
and (c) show the time dependence of the fluctuation for Ti(t) 
and £,i(t), respectively. The heat current into the island cor- 
responding to this fluctuation is plotted in (d) , and (e) shows 
the charge current flowing through the island. 



the curve quickly depart from equilibrium to either higher 
or lower temperatures. The solution of the saddle-point 
equations follows A before the measurement and B after 
the measurement (Figs. [5]b-e). 

Since S c = 0, the only contribution to path integral 
([3]) comes from the island term C£/T 2 and is evaluated 
as 



"tt 2 






2tt 2 y T « 


.3(5/ j 


f £_T 2 dt 


= exp 


3(5/ J Ta 



P(T e )=exp — / iT 2 dt =exp -— / T£f(T)dT 

(9)' 

Thus, in order to find P(T e ), we only need a function 
£f (T/) satisfying S c (£f (!», T/) = at branch A. 

The connector action can generally be written in the 



form 

x (e-xr-fr* _ i) + _ /^(eXJ+fec „ 

(10) 

with f a/I = {cxp[(e - fia/i)/T a /j] + l} -1 . 

To prove the validity of the method for the equilibrium 
case, let us set all the chemical potentials to and all the 
reservoir temperatures to T a . This implies fir = xi = 0. 
Using the fact that for a Fermi function / = — e e//T (l — /), 
we observe that S c = regardless of contact properties 
provided £/ = £f (I» = I/Tl-1/Tj. Substituting this to 
Eq. ([9]) reproduces the equilibrium distribution, Eq. fTJ. 

Out of equlibrium, the further analytical progress can 
be made in the case when the connectors are ballistic, 
T n = 1. Such a situation can be realized in a chaotic 
cavity connected to terminals via open quantum point 
contacts. The connector action reads 91, 



2fi aX i + T aX j + [vr 2 T 2 /3 + fi\ £/ 



1 - T c ^i 



1 + T/Cr 



(11) 



Let us first assume two reservoirs with Ti = T% = Tt,. In 
this case the general saddle-point solution for the poten- 
tial follows from Kirchoff law: /!/ = (g/J,i + A*2)/(l + g) 
with g = Gl/Gr. For the charge counting field we get 
Xi = — The most probable temperature T a is given 
by T 2 = T 2 + 3 5 ( M i - M2) 2 /k 2 (l + ff) 2 L and function 
£f(Ti) is expressed as 
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(12) 



Ti(Tl,Ti + T%) 
Substituting this to Eq. (J9j> yields for the full probability 
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(13) 



In the strong nonequilibrium limit V = (/ii — ^2) 3> Pl, 
i.e., T a » Tf, this reduces to 



Pbaii oc exp 



27r 2 fc B (T e + 2T a )(T e -T a ) 5 



3(5/ 



3T 2 



(14) 



The logarithm of this probability is plotted as the lower- 
most line in Fig. [3l 

If the island is biased by temperature difference, T\ = 
Tl >• T 2 , V = 0, the probability obeys the same Eq. ff 
withT 2 = ff T 2 /(l + . 9 ). 
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FIG. 3: (Color online.) Logarithm of TFS probability 
P(T e ) in a few example cases. Solid lines from top to 
bottom: temperature bias with symmetric tunneling con- 
tacts, T a = Ti/%/2, T% = (magenta); Gaussian equilib- 
rium fluctuations (black), nonequilibrium fluctuations with 
T a = V3\eV\/(2irk B ), Ti = T 2 = for symmetric tunneling 
and ballistic contacts (blue and red lines, respectively). The 
dashed lines are Gaussian fits to small fluctuations (T e —T a ) <C 
T a , described by the heat current noise Sq at T e ~ T a . 



For general contacts, the connector action and its 
saddle-point trajectories have to be calculated numeri- 
cally. For tunnel contacts, the full probability distribu- 
tion is plotted in two regimes in Fig. [3] The distribution 
takes values between the ballistic and equilibrium cases. 
Let us understand this concentrating on Gaussian regime 
and inspecting the variance of the temperature fluctua- 
tions for various contacts. This variance is related to the 
zero-frequency heat current noise Sq via 



2G th C(ST 2 ) =S Q = dp c \^ 



(15) 



In equilibrium, S 



(eg) 



2G th T 2 by virtue of the 



fluctuation-dissipation theorem. For an island with equal 
ballistic contacts driven far from equilibrium, V ^> Te,, 
S^ 1 = V3GV 3 /(8tt) = Gth(T„)2* i.e., only half of S Q eq) . 
The reduction manifests vanishing temperature of the 
reservoirs. Most generally, for contacts of any nature, 
the heat current noise reads 



V4 £9) = 2 



1 



a Q ^ Fa i 



(16) 



where F a = £„T„ Q (1 - T n «)/Ej„" is the Fano factor 
for a contact a, aQ w 0.112 being a numerical factor. For 

two tunnel contacts we hence obtain S^ n ~ 0.723S^ q \ 
a value between the ballistic and equilibrium values. For 
contacts of any type, the variation of temperature fluc- 
tuations is between the ballistic and tunneling values. 



For rare fluctuations of temperature, \T e — T a \ ~ T a , 
the probability distribution is essentially non-Gaussian 
in contrast to the equilibrium case. The skewness of the 
distribution is negative in the case of voltage driving: 
low-tempcraturc fluctuations (T e < T a ) are preferred to 
the high-temperature ones (T e > T a ). In contrast, bi- 
asing with a temperature difference (uppermost curve in 
Fig. [3} favours high-temperature fluctuations. 

The non-Gaussian features of the temperature fluctua- 
tions can be accessed at best in islands with a large level 
spacing, that is smaller than the average temperature say, 
by an order of magnitude. Many-electron quantum dots 
with spacing up to 0.1 K/fc^ seem natural candidates for 
the measurement of the phenomenon. The most natural 
way to detect the rare fluctuations is through a threshold 
detector 12( , which produces a response only for temper- 
atures exceeding or going under a certain threshold value. 
Besides the direct measurement of temperature, one can 
use the correlation of fluctuations. For example, Fig.f^e) 
shows that the fluctuation of the temperature also causes 
a fluctuation in the charge current. Observing the latter 
may thus yield information about the former. 

To conclude, we have evaluated non-equilibrium tem- 
perature fluctuations of an example system beyond the 
Gaussian regime. The method makes use of saddle-point 
trajectories and allows to describe electric contacts of ar- 
bitrary transparency. 
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